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Abstract 

We show that the colour-singlet process T — > J/'^ + ccg significantly 
contributes to the inclusive process T J/^ + X. We calculate the par- 
tial width and the momentum distribution of the produced J/^ in this 
channel. The obtained width is comparable to the experimental data. The 
momentum distribution is fairly soft, which is in contrast to the results 
from the colour-octet processes discussed earlier in literature but is con- 
sistent with the CLEO measurements. Further experiments to check this 
contribution are suggested. 



*e-mail: tpc@sdu.edu.cn 



1 



Heavy quarkonium decay and production are traditional arena for QCD. The 
descriptions |T| of such processes are nearly standard as a factorization proce- 
dure: The annihilation of bb or cc to gluon(s) or its creation from gluon(s) can 
be calculated via perturbative QCD, since the heavy quark mass sets a hard 
scale for corresponding processes. The bound state of the heavy quarkonium 
contains non-perturbative effects. This part is 'factorized out' and is described 
by QCD-inspired models 0. Such a procedure is successful in describing many 
experimental data. However, recently a large discrepancy for the prompt produc- 
tion of J/\I' and has been observed in Fermilab Tevatron between the data 
and the colour-singlet model prediction. This lead to a re-inspection of the de- 
scription of heavy quarkonium decay and production (for a review, see |^). It was 
found that the 'factorization' in colour- singlet model had never been proved. And 
the colour-octet mechanism which is based on the NRQCD factorization scheme 
(see, e.g. [0) was proposed as a complement to the colour- singlet mechanism to 
explain the above mentioned discrepancy. It has been shown that the colour-octet 
mechanism plays an important role in many processes. For example, it has been 
argued that two processes T ^ J/"^ + gg ^ (see Fig. 1(a)) and T — > J + g 
[0 (see Fig. 1(b)) where the colour-octet mechanism plays the key role provide 
the leading contribution to the inclusive process T J/'^ + X. Experiments by 
the CLEO collaboration reported a branching ratio of (1.1 ± 0.4 ± 0.2) x 10^^ 
and quite a soft J/\l/ momentum spectrum (see dots in Fig.3) P|. The calculated 
branching ratios of these two colour-octet processes are indeed comparable to 
the data. However, the J/\l/ momentum spectra of both processes peak at the 
largest value ~ 4.2Ge\^/c, which is in strong contradiction to the CLEO data. 
(For T J/^ + gg, see the dotted line in Fig.3; for T J/^ + g, the J/'^ 
momentum spectrum is a ^-function with the pole at ~ A.2GeV/c.) 
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In this note, we show that a colour- singlet process, T J/\E' + ccg (see 
Fig. 2), also significantly contributes to the inclusive J/\l/ production in T decay, 
and that the J/"^ momentum spectrum from this channel has to be soft, which is 
consistent with the CLEO data. The latter point can be seen from the following 
simple kinematic analysis. We recall that angular momentum conservation and 
charge conjugation invariance demand that T decay mainly through three (real 
or virtual) gluons. It is clear that T can decay via one real and two virtual gluons, 
and both of the virtual gluons split into colour-octet cc (see Fig. 2). In this case, 
a c and a c originated from different gluons can be in the colour-singlet, and they 
can form a J/\E' when they are in the proper angular momentum state and their 
invariant mass is near J/\l/ mass shell. Generally, the remaining cc created in 
this process hadronizes into open charm particles, so the final particle system 
excluding J/\l/ has a minimum invariant mass 2M£,o. This implies that the J/\l/ 
momentum has a maximum of 3.3GeV/c in this channel. From the Feynman 
diagram, we see also no indication that the J/\l/ momentum spectrum peaks at 
this largest value, so it should be soft. The exact shape depends on the details 
of the dynamics which is discussed later in this paper. We also note that if only 
one instead of two cc is created in the 3-gluon decay channel, which to the lowest 
orders in as correspond to the above mentioned two colour-octet processes p, , 
the J/\l/ momentum spectrum is very hard and in strong contrast to the CLEO 
data, as cited at the end of last paragraph. Hence we conclude that in the 3-gluon 
decay channel and to the lowest orders in a^, the process in Fig. 2 is the only case 
where the J/"^ momentum spectrum is soft. 

The subsequent question is then: how large is the branching ratio of this 
colour-singlet process? At first sight, one may expect that it is much smaller 
than that of the colour-octet process T J/ip + gg discussed in 0, since in 
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the perturbative (short-distance) phase the colour- singlet process has one more 
as and one more virtual gluon propagator, both of which give more suppression. 
But a careful inspection tells us that it can be of the same order of magnitude 
as that from the colour-octet process T — > J/ip + gg, thus comparable to the 
data 1^. The reason is as follows. It is true that in the perturbative phase 
the colour-singlet process we consider receives more suppression than the colour- 
octet process 0. But in the non-perturbative (long-distance) phase, where the 
colour-singlet/octet cc form the J/ip resonance with the probability described 
by the colour-singlet/octet matrix elements [^, the much larger colour-singlet 
matrix element provides an enhancement. Using the value of 0, the colour-singlet 
matrix element < 0'(^'^(^Si) > is about 52 times as large as the colour-octet one 

we use the more updated value in the ratio — j/^ 
can even be 210 to 360. This enhancement can compensate the perturbative 
suppression. So we expect a non-negligible partial width of the colour-singlet 
process T cc{^Si, 1) + ccg J/\E' + ccg, where we use {^Si, 1) to denote the 
quantum numbers of angular momentum and colour for cc. Furthermore, if the 
invariant mass of the colour-singlet pair cc(^S'i) is near the mass shell of ^'(25'), it 
can also form a \E'(2S') resonance, which decays into J/\I' with a branching ratio 
of 54.2% |T0[. This indirect contribution T \E'(2S') + ccg J/\E' + X provides 
another source of the inclusive J/'^ production in T decay. 

Encouraged by these qualitative analyses we calculate the partial width of 
T and momentum distribution of J/\E'(\I'(2S')) for the direct process shown in 
Fig. 2 and the process T cc(^Si, 1) + ccg \E'(2S') + ccg. Since both of them 
are colour- singlet processes and the charm quark is heavy, the treatment for the 
bound states can be a conventional wave-function approach where the relative 
momentum between c and c is vanishing, namely same as the case of positronium 
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see, e.g. [|rT|). Similar to |]TT|, the differential width of the direct process (see 



Fig. 2) can be formulated as 

dT _ \Bj/^Br < cc{^Si,l)ccg\S\bb{^Si,l) > ^ 



(1) 



dR T 
where dR is the 8-dimensional phase space volume element for J/\l/ and c, c,g; 
S is the S-Matrix; Bj/^ and i?x are related to the origin values of the wave 
functions of J/\l' and T: 

^ (2) 



Br = (3) 

For convenience, we normalize all particle states to 2EV (where E is the particle's 
energy and V is the volume of the total space) except for T, which is normalized 
to 1. This is why the coefficients of the wave functions in (2) and (3) are different. 
In (1) the sum over all spin states for final particles and average of the 3 spin 
states for T are not explicitly shown and the 'time' T is 27r5(0). The corresponding 
expression for the differential width of the process T cc(^Si, l)+ccg — > \I'(2S') + 
ccg is the same. The only differences lie in the phase space and the wave function 
at the origin, where the corresponding quantities for \E'(2S') should be used to 
replace those of J/\&. 

For the perturbative sub-process bb{^Si, 1) g*g*g —>■ cc cc g, we consider, as 
usual, only the six lowest-order Feynman diagrams at the tree- level (see Fig. 2). 
The treatment of the infrared behaviour of the real gluons in these diagrams is 
simple. It is much similar to that in the T —>■ 3g process (see, e.g. [H). For the two 
diagrams where the real gluon vertex lies between the two virtual gluons, there 
is no infrared singularity due to 'controlling momentum' [Q. For the rest four 
diagrams, the infrared gluon in each diagram leads to an eikonal vertex and an 
eikonal quark propagator [^]. These four diagrams form two pairs. In each pair. 



the only difference between the two diagrams is the position of the real gluon ver- 
tex. The eikonal vertices in them have opposite signs, so the infrared singularities 
of the two diagrams cancel with each other. Summing over all the six diagrams, 
we obtain the infrared safe matrix element < cc{^Si, l)ccg\S\bb{^Si, 1) >, which 
is needed in calculating ^ from Equation(l). 

To get the numerical results, we need to know the value of a.^ and those of 
the \E'j/^(0) and \E't(0) (see Equation(2),(3)). In order to compare the results 
with those of the colour-octet processes [|, 0], we use the same values for these 
parameters. Especially, the scale of is Mj/q,. The square of the radial wave 
function at the origin for "${23) have been calculated using different potentials 
[0, Consistent with these calculations, we use [^*|^g^(o)}2 ~ 1/2 approximately, 
where TZ represents the radial wave function. The quark masses are taken as 



rrir 



M 



and TTih = The partial width of the direct production process 



T cc{'^Si, 1) + ccg J/^ + ccg is calculated to be 28.4eV. Normalized by 



the full width of T |jTO[, we obtain the corresponding branching ratio of about 
0.54 X 10"'^. The indirect contribution to the branching ratio from \E'(2S') is small 
~ 0.05 X 10^^. So the branching ratio from these two sources is 0.59 x 10~^, 
which is consistent with the CLEO result within the error. 

If we take sum of the branching ratio obtained above for the colour-singlet 
process we consider, and those for the colour-octet processes discussed in 0, 
we obtain the total branching ratio of 1.21 x 10^'^. This result is still in agreement 
with CLEO data (but is larger than the ARGUS upper limit). Because the value 
of the colour- singlet matrix elements can be extracted from the pure leptonic 
decays of T and J/\E' etc., which are electromagnetic process and where only 
the colour-singlet mechanism works, the uncertainty is rather small. This is in 
contrast to the case of extracting the value of the colour-octet matrix elements. 
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The latter can only be extracted from J/\E' production experiments. It is therefore 
very sensitive to whether all the colour- singlet processes have been considered. 
Furthermore, the perturbative part is sensitive to the value of a^. It is unclear 
whether as{Mji^,) or as{M-^) should be used 0. If ^^(Mt) instead of ^^(Mj/^) 
is used in both the colour-singlet process we considered here and the colour-octet 
processes considered in 0], the branching ratio of the colour-singlet process 
decreases to 0.12 x 10~^ which is still of the same order of magnitude as that of 
the colour-octet processes. But the theoretical total width is smaller than the 
CLEO data. It is now impossible to judge which mechanism is more important 
if we only look at their contributions to the branching ratio. 

In contrast to the branching ratio, the shape of the calculated J /'^ momentum 
spectrum does not depend on the values of as and the non-perturbative matrix 
elements. It is not influenced by the uncertainties coming from these parameters. 
Hence it can act as a 'label' for a given decay channel. The J/\E' momentum 
distribution for the direct colour- singlet process(see Fig. 2) is shown in Fig. 3 by 
the solid line. It has a peak near 2.2GeV/c. The shape is similar to the CLEO 
data. About 90% of the total events reside in the range smaller than 3 GeV/c, 
which is also consistent with the CLEO data. The momentum spectrum of \I'(2S') 
is shown by the dashed line, and the magnitudes of the differential width have 
been multiplied by the branching ratio of \E'(2S') ^ J/\E' + X' . In the process 
*(25') J/^+X', since Mvi>(25)-Mj/vi, << Mj/^, J/^! carries most of the energy 
of \E'(2S'). Hence the shape of the corresponding J/\E' momentum spectrum is 
approximately the same as the '^{2S). We see from Fig. 3 that the spectrum shape 
in the indirect process is similar to that in the direct process. For comparison, 
we also show in Fig. 3 the momentum spectrum of the colour-octet process (see 
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Fig. 1(a)) by the dotted line Q. We see clearly that the CLEO data favour the 
colour-singlet process considered in this note. 

Finally, we would like to emphasize that it is not difficult to check experimen- 
tally whether the colour-singlet process we consider indeed plays an important 
role for T J/'^ + X. This is because the J/\E' in this colour- singlet process 
is definitely associated with open charm particles and its momentum is less than 
3.3GeV/c. 

In summary, both qualitative analyses and quantitative calculation show that 
the colour-singlet process T J/"^ + ccg gives the branching ratio comparable to 
experimental data and a similar soft J/\l' momentum spectrum as measured by 
CLEO Collaboration |Q as well. Among the processes which have been studied 
theoretically up to now, this is the only one which has such a feature. This seems 
to suggest that this process provide the main contribution to the inclusive J/"^ 
production in T decay. 
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^We use Equation (24) in ^ to draw the line. We have changed the variable to J/^* 
momentum and normahzed the expression on the right hand side of the equation to the width 
for this channel 
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Figure Captions 

Fig.l Diagrams illustrating the colour-octet processes T ^ J/\& + (7(7 @ (a) 
and J/^ + ^ (b). 

Fig. 2 One of the six diagrams for the direct colour-singlet process T — >■ 
cc{^Si, 1) + ccg — > J/\E' + ccg. The real gluon vertex can be between those of 
the two virtual gluons or in either side of them. And the two virtual gluons can 
change their order. 

Fig. 3 The momentum spectrum of J/\E' in the direct colour-singlet process 
we consider (solid line) and the momentum spectrum of ^(25') multiplied by the 
branching ratio of \E'(25') J/'^ + X' for the process T — > \E'(25') + ccg (dashed 
line). The CLEO data p| is shown by dots. The result of the colour-octet process 
of Fig. 1(a) is shown for comparison (dotted line). 
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